Atom-Molecule Laser Fed by Stimulated Three-Body Recombination 
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Using three-body recombination as the underlying process, we propose a method of coherently 
driving an atomic Bose-Einstein condensate (BEC) into a molecular BEC. Superradiant-like stim- 
ulation favors atom-to-molecule transitions when two atomic BECs collide at a resonant kinetic 
energy, the result being two molecular BEC clouds moving with well defined velocities. Potential 
applications include the construction of a molecule laser. 

PACS numbers: 03.75.Mn,03.75.Pp,34.10.+x 
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One of the major goals of the recent research in BEC 
physics is the creation of a molecular condensate 0. 
Many applications are envisioned for such a system, in- 
cluding coherent quantum chemistry 0. One approach 
to forming a molecular BEC utilizes an atomic conden- 
sate along with a method of coherently driving pairs of 
atoms into (bound) molecular states. Two mechanisms 
have been proposed for this: (i) photoassociation, in 
which an external laser couples atomic states to molec- 
ular states [3, and (ii) control of the scattering length 
via external magnetic field ramps tuned near a Feshbach 
resonance ^ . Another problem of interest to the BEC 
community is the realization of an atom laser Q . A pro- 
posed technique for this is the application of a time- 
dependent light pulse resulting in Bragg scattering of the 
condensate cloud, in which a fraction of the condensate 
receives a large momentum kick and attains velocity rel- 
ative to the parent cloud. 

In this Letter we propose a new mechanism for the for- 
mation of a molecular BEC from an atomic BEC. This 
mechanism relies on a twofold bosonic stimulation of the 
three-body recombination process — both by the number 
of atoms present and by the number of molecules. We 
show that two atomic condensates with a well-defined 
(resonant) relative velocity will mutually stimulate tran- 
sitions into molecular states, resulting in the production 
of two coherent molecular clouds, each with a well-defined 
momentum. 

Consider first the elementary process involving only 
three isolated atoms. The atoms, initially in an unbound 
state, are converted through recombination into a bound 
molecular state with 1/3 of the recoil energy, while a 
third atom carries away the remaining 2/3. Conserva- 
tion of momentum requires that the resulting atom and 
molecule have momenta of equal magnitude (hko), but 
in opposite directions. The value of fco is fixed by con- 
servation of energy to be fcg = y^AmEo/ifP, where Eq is 
the molecular binding energy and m is the atomic mass. 
This process is characterized by a transition matrix ele- 
ment Tif between the initial and final three-body states. 
The spontaneous decay of an ultracold atomic sample 
of density n through three-body recombination is usu- 
ally characterized by a loss rate coefficient K^, such that 



dn/dt — —SK^n^ /6. Here is the number of atomic 
triplets available for recombination in the considered vol- 
ume and the factor of 3 denotes the atoms lost in each 
recombination. An additional factor of 1 /6 multiplies the 
right-hand side if all three atoms share the same state, 
i.e., if we consider the decay of a pure condensate 0. 
The connection between Ti f and K3 is given by 
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If the recombination forms a weakly bound level with 
Eq « /{raa^) then fcp = 2/(-\/3a), where a is the s-wave 
scattering length. In this case, working in the ultracold 
limit and using momentum normalized states, Ref. [8] ob- 
tains Tif — C,h^ a^/'^ /m 0, where C is a numerical con- 
stant. More generally has been shown to scale with 
the fourth power of the scattering length in most cases, 
which for the weakly bound level 
gives (3 = (12C^)/(7r\/3). Ref. % finds /? « 23, while a 
more rigorous calculation, valid for a wide range of in- 
teratomic potentials Q, yields f3 ~ 180 for a > and 
P « 1014 for a < 0. (We omit a sin^ modulation factor 
discussed in [9,10].) 

We are interested in the process in which three-body 
recombination is stimulated by bosonic enhancement. 
We consider two atomic condensate clouds moving with 
relative velocity hko/m, and choose a reference frame in 
which the two clouds move with equal velocities towards 
each other (see Fig. ^ along the direction k. Let the 
wavevector of the atoms moving to the right (left) be 
denoted by ki (k2) and the population of the mode by 
ni (712). Note that ki = — k2 = ko/2k. Recombination 
of three atoms with momentum Kki can lead to a final 
state in which the resulting atom has momentum h'k2, 
in which case the resulting molecule has the momentum 
Kkyni — 2hko k. This transition will be strongly favored 
due to bosonic stimulation if the mode k2 is highly pop- 
ulated. We recall the fact that the magnitude of the mo- 
mentum is fixed by energy and momentum conservation 
and only its direction has an arbitrary value. Therefore, 
it is possible that the atom-molecule pair is ejected along 
a direction different from k. However, this process does 
not benefit from bosonic stimulation. By symmetry, re- 



combination of three atoms with momentum ?ik2 will be 
strongly enhanced if the resulting atom-molecule pair is 
ejected along the direction k. In this case the momentum 
of the molecule has the value ^^^2 — —2 hk^ k while the 
resulting atom has momentum fi,ki . We denote the num- 
bers of molecules with wavevectors kmi and ki„2 by n^i 
and nm2, respectively. Finally, note that, for any three- 
body recombination the reverse process is also possible. 
For example, an atom with momentum h\i2 that collides 
with a molecule with momentum /ikmi can induce it to 
break up, resulting in three atoms with momentum ftki . 
In summary, the four modes considered here are coupled 
by the transitions 

3A(ki) ^ A(k2) + A2(k^i), 

3A(k2) ^ A(ki) + A2(k^2). ^' 

The many-body interaction Hamiltonian describing 
the reactions in Eqs. is 

Hi = huj (^cl c]^;^ Ci Ci Ci + c\ 4^2 C2 C2 C2) + h.c. (2) 

Here ci_2 (cmi,m2) are the annihilation operators for the 
two atomic (molecular) modes and the coupling fre- 
quency is given hy huj = Tif/V'^^'^, where V is the 
quantization volume. For a weakly bound molecule 
w = [Ch/{ma'^)]{a^/V)^/'^. The two operators, 

N ^ c\ci + C2C2 + 2 cjjji Cml + 2 c]jj2 Cm2 , 

= Ci - C2 -I- 4 cl^^ Cml - 4 cl,2 Cm2 , 

commute with the Hamiltonian in Eq. (0) as a conse- 
quence of the conservation of the total number of atoms 
and of the total momentum. 

We now derive the equations of motion for the creation 
and annihilation operators in the Heisenberg representa- 
tion. Adopting the notation ai^2iT) and 61.2 (t) for the 
atomic and molecular Heisenberg operators divided by 
■y/iV (in which fast oscillations with the frequency asso- 
ciated to the non-interacting energy were factored out), 
the equations of motion read. 
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We use a scaled time t = Vlt with O = N'^/"^ uj — r?!"^ Ti j 
[n = N/V is the total concentration of atoms). For a 
shallow molecular level = [C,h/ {ima?)]{na?Y/'^ . Equa- 
tions describe the dynamic population transfer 




FIG. 1: Schematic diagram of the four possible states of 
the atoms: atomic states with momenta ftki and h\i.2 have 
populations rii and n2 ; molecular (molecular mass 2m) states 
with momenta fikmi and ftkm2 have populations nmi and nm2 • 



between the atomic and molecular modes due to reso- 
nant three-body recombination. They are based on the 
assumption that all two-body elastic collisions (atom- 
atom, molecule-molecule and atom-molecule) and also 
coUisional loss of the weakly bound molecules to deeper 
molecular states can be neglected. In addition, we have 
neglected the spatial and temporal changes of the slowly 
varying envelopes of the condensate clouds. The sim- 
plest approximation to the above equation consists of 
replacing the operators by complex numbers equal to 
their averages: ai,2(0 = (01,2 (0) and 61,2 (*) = (&i,2(0): 
i.e., the classical field approximation. The resulting sys- 
tem of ordinary differential equations, with initial con- 
ditions, can be integrated numerically. We consider the 
experimentally relevant case of an initial state with zero 
molecules. We also assume that initially the N atoms 
are evenly split between the two momentum eigenstates: 
A^i.2 = ni^2{t = 0) = A^/2. Assuming real initial con- 
ditions, we set 61,2(0) = and 01.2(0) = l/\/2- The 
approximate version of the system of Eqs. can 
be solved to describe the population transfer from the 
atomic to the molecular modes in a process characterized 
by the time constant Tr = 2tt/VI. A characteristic that is 
similar to the superradiance phenomenon is double am- 
plification (stimulation): after an initial buildup period 
for the molecular population, rates of further transitions 
are proportional both to the number of existing molecules 
in the recombined momentum state and to the number 
of atoms in their relevant momentum state. 

The question of the validity of using the classical field 
approximation can be studied in general by accounting 
for higher-order terms in a cluster expansion of the four- 
operator products in Eqs. ®-® (see, e.g., Ref. 11]). 
However it is sufficient for our purpose to compare the 
results of the approximate version of Eqs. lO^® to 
the many-body Schrodinger equation written using a ba- 
sis of Fock states. We designate these basis states by 
|ni, ?i2, 71ml, nm2) where ni,2 (n-mi, 111,2) are the occupa- 
tion numbers of the two atomic (molecular) modes. Due 
to the two conservation laws ni = iVi -f n^2 — 3nmi and 
n2 ~ N2 + n-aii — 37T,m2, wc Can restrict the basis set to 
only those accessible by three-body recombinations from 
our initial state and thus use only two occupation num- 
bers, say riiiii and nm2, to parameterize the Fock states. 
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FIG. 2: Comparison of the scaled-time-dependent average 
number of atoms (solid lines) and molecules (dotted lines), 
obtained in two ways: by solution of the approximate Heisen- 
berg equations (thick lines) and by solution of the many-body 
Schrodinger equation (thin lines). The three panels corre- 
spond to different total numbers of atoms A^: top — 64, 
middle N = 640, bottom = 6400. 

With this notation, the many-body wavefunction is 

>• (7) 

"ml ,nm2 

After eliminating the fast osciUations and introducing the 
time scaling t' = ojt, the Schrodinger equation becomes 

= V ('^1 + + 2)(ni + 3)n2nml /n„i-l,n„2 

+ \/ni{ni - l)(ni - 2){n2 + l){n^i -f 1) /„„i-i-i,«,„2 

+ ^("2 + l)(n2 + 2)(n2 -I- 3)ninin2 /n„,i,n„2-l 

+ sjniini - l)(n2 - 2)(ni + l)(n„i2 + 1) /n„,i,«„2+i- 

Solving these equations with the initial conditions 
/o,o(0) = 1 and /n„i,n„2(0) = for n,„i ^ and 
?^m2 7^ 0, allows us to calculate the average numbers 
of atoms populating any of the four considered modes. 
These averages are shown in Fig. |21 for different values 
of Ni — N2 = N/2. Comparison with the corresponding 
results obtained using Eqs. shows agreement over 

the initial period of fast transitions of atoms into molecu- 
lar states. As expected, the classical field approximation 
fails at shorter times for lower values of N, while for 
sufficiently large N it can be used to accurately describe 
the initial transition burst converting atoms to molecules. 
Equation © allows us to describe the population transfer 
within the four-mode approximation without making any 
additional approximations. However, calculations using 
the Fock basis become prohibitive if more than a few 
thousand atoms are considered. 

Experimental arrangements that use the process de- 
scribed above to produce coherent molecular clouds can 



FIG. 3: Solutions of Heisenberg equations with elastic loss 
(top panel: 7 = 0.1, bottom panel 7 = 1.5). Solid line: 
atomic population fraction |5i^2(''")P, dotted line: molecular 
population fraction |bi,2(T")p. 



be readily envisioned. The simplest experiment would 
begin with a time-dependent Bragg splitting of a single 
atomic cloud as in the experiment of Kozuma et al. (||, 
but with the particular choice of momentum kick that 
matches the three-body resonance condition. In addition, 
the value of the molecular binding energy can be tuned 
using an external magnetic field [l3.fl3| . In order to take 
full advantage of the stimulated transition into molecular 
states, (i.e., to transform a large fraction of the atomic 
clouds into molecular clouds) the collision time of the two 
atomic clouds should be comparable to the characteristic 
time of the three-body recombination process icoii ~ Tfl- 
Here icon ~ Rc/if^-ko/m) and Rc is the characteristic size 
of the condensate along the collision direction. Consider 
a sample of ^^Rb with a concentration of lO^^/cm"^. As- 
suming that the scattering length is adjusted with the 
help of a Feshbach resonance to a = 400 a. u. and that 
Eq — ft^/ma^ then we estimate Tr = 0.45 ms (using 
K3 given by j^). The resonant process described here 
is characterized by a finite width if additional processes 
that limit the lifetime of the involved states are consid- 
ered. One such process is the coUsionally induced decay 
of the molecular level to lower molecular states. Only 
a very rough estimation of the lifetime is available for 
molecules formed by atomic species commonly used in 
BEC experiments. Some estimations (see, e.g, cite 
values of at least 1 ms for this lifetime at molecular den- 
sities of 10^^/cm'^. In addition, recent experiments in 
^^Rb suggest that this is a conservative estimation, 
and that the actual lifetime might be somehow longer. 

Up to this point we have neglected elastic collisions be- 
tween atoms and molecules. In fact, if the resonant value 
of the relative velocity is not enforced by the experimen- 
tal setup, depletion of the atomic clouds by elastic atom- 
atom scattering is the main loss process ta king place be- 
fore the two condensate clouds separate [ij- As dis- 
cussed in Ref. 0, elastic scattering must be accounted 
for beyond the mean-field (i.e., Gross-Pitaevskii) approx- 
imation in order to observe the condensate depletion due 
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to elastic scattering. We can include this process in our 
simple model as a loss term in the interaction Hamil- 
tonian. Considering only atom-atom collisions involving 
one atom with wavevector ki and one with k2, Eqs. Q 
and I0J in the classical field approximation become 

i—ai — Salal bi 0.2 + ^2 ^2 0,2 02 — «7 |a2|^ ai , 

1—0.2 — ioXaX 62 oi + b\ hi ai ax ~ 27 laiP 02 . 
ar 

Here the parameter 7 is the elastic scattering rate 
constant divided and is given by 7 = 

{hko / m) an / {2hil) 15], where a is the elastic scattering 
cross section at the relevant two-body collision energy 
^'coii = h^kQ/m. Eqs. © show that the two processes — 
resonant three-body recombination and condensate loss 
due to elastic scattering — compete with each other and 
that the value of the parameter 7 is crucial in determin- 
ing which process dominates. Fig. O shows the solution 
of the Heisenberg equations that include the elastic loss 
terms for two extreme cases. The example with 7 = 0.1 
illustrates the case of little elastic loss allowing almost 
70% of the atoms to be transformed into molecules if 
icoii ~ Tfi. However, if 7 = 1.5, for the same collision 
time more than 80% of the atoms are lost due to elas- 
tic collisions and only about 8% end up in the molecular 
clouds. 

The simplest estimation of 7 can be attempted in the 
low energy limit (fco ^ where a ~ 87ra^. This 

leads to 7 = (87r/-\/3C)(na'^)~^/^ and yields high values 
of 7 even for very large values of the diluteness parameter 
(i.e., na^ = 1 corresponds to 7 = 1.5). However, note the 
importance of using an exact value of a corresponding to 
the correct Eco\i- For example, in the case of a weakly 
bound molecular level [£^0 ~ fi? / {ma?)], ko « 1/a, which 
is beyond the validity of the low energy approximation. 
Apparently, changing the value of the parameter na^ by 
changing either n or a is the primary way of shifting 
the balance between the two processes. However, spe- 
cial conditions might enhance or suppress either process. 
One such situation may occur when a approaches zero 
and therefore elastic scattering vanishes while three-body 
recombination still has a considerable probability. Cal- 
culations by Esry et al. predict a significant K3 for 
a — 0. Similarly, Stenger et al. Q| measured three-body 
recombination loss from Na condensates in the vicinity 
of Feshbach resonances, observing a large amount of loss 
even at small values of a. Consider the Feshbach reso- 
nance located near B = 535.7 G in ^^Rb. Tuning the 
magnetic field such that a = {B = 537.6 G) leads to 
Eq = 75 ^K, and for a density of lO^^/cm'^ we obtain 
K3 = 5.4 X 10~'^^cm^/s (as estimated in Ref. 0), giving 
Tr « 28 ms. For Na with a = (B = 907G), Ref. 
measures K3 w lO^^^cm^/s, giving Tr w 0.24ms. Al- 
ternatively, one can choose the experimental conditions 
such that the elastic collisions are strongly suppressed 



exactly at the resonant value of the collision energy i^coU 
due to the presence of a Ramsauer minimum (see, for 
example, [l6j). The exact resonant matching between 
the position of the Ramsauer minimum and that of the 
binding energy of the molecular state (i.e., i?coii = Eq/Z) 
can be achieved near some Feshbach resonances. For the 
above mentioned resonance in ^^Rb, this matching occurs 
for B = 547 G, when a = -353 a.u., and Eq = 1.05 mK. 
In this case, considering n = lO^^/cm"^, one finds Tr « 
1.3 ms {K3 = 9.2 X lO-^^cm^/s cf. 0). Since producing 
an atomic condensate at negative a might be a problem, 
we assume here that it is produced at some positive a and 
that the magnetic field is swept into resonance precisely 
when the two atomic clouds start colliding. 

In conclusion, resonant three-body recombination that 
benefits from the superradiant-like bosonic stimulation 
could be a promising candidate for producing coherent 
samples of molecules from existing atomic condensates. 
The simplest experimental setup that takes advantage 
of this process involves a collision of two atomic con- 
densates, which produces two molecular condensates in 
counterpropagating momentum eigenstates. This setup 
could be particularly promising for the construction of a 
molecule laser by colliding two atomic condensates. 
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